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ABSTRACT 



Diagnostics of electron temperature (T e ), electron density (n e ), and hydro- 
gen column density (Nh) from the Chandra High Energy Transmission Grating 
spectrum of He-like Ne IX in TW Hydrae (TW Hya), in conjunction with a clas- 
sical accretion model, allow us to infer the accretion rate onto the star directly 
from measurements of the accreting material. The new method introduces the 
use of the absorption of Ne IX lines as a measure of the column density of the 
intervening, accreting material. On average, the derived mass accretion rate for 
TW Hya is 1.5 x 10~ 9 M Q yr _1 , for a stellar magnetic field strength of 600 Gauss 
and a filling factor of 3.5%. Three individual Chandra exposures show statisti- 
cally significant differences in the Ne IX line ratios, indicating changes in N H , 
T e , and n e by factors of 0.28, 1.6, and 1.3, respectively. In exposures separated 
by 2.7 days, the observations reported here suggest a five- fold reduction in the 
accretion rate. This powerful new technique promises to substantially improve 
our understanding of the accretion process in young stars. 

Subject headings: accretion, accretion disks — stars: formation — stars: individual 
(TW Hydrae) — techniques: spectroscopic — X-rays: stars 
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1. Introduction 

In the standard scenario for magnetic accretion onto classical T Tauri stars (CTTS), 
material accelerates along magnetic field lines from a truncated accretion disk toward 
the star, producing an X-ray emitting shock of a few MK near the stellar surface. Many 
accretion-related phenomena observed in ultraviolet and optical spectra, e.g. excess 
ultraviolet continuum and the filling in or "veiling" of photospheric absorption lines by hot 
continuum, are interpreted as the stellar atmosphere's response to the X-ray shock. Optical, 
near infrared, and ultraviolet emission lines also show accretion-related broadening of a few 
hundred km s" 1 . Optical spectral signatures facilitate a number of indirect techniques for 
estimating the mass accretion rate and hot spot filling factor. High resolution X-ray spectra 
of CTTS with the Chandra and XMM-Newton gratings provide line ratio diagnostics of the 
electron temperature (T e ) and electron density (n e ) in the shock which supports this basic 
magnetospheric accretion model. Diagnostics of the X-ray shock in principle can provide a 
more direct measure of the accretion parameters (e.g. Giinther 2011), but they have often 
significantly underestimated the accretion rate relative to more traditional methods (e.g. 
Calvet & Gullbring 1998; Hartigan & Kenyon 2003). 

Kastner et al. (2002) presented the first evidence for high n e in a CTTS using the 
density-sensitive He-like Ne IX diagnostic ratios from a Chandra High Energy Transmission 
Grating (HETG) spectrum of TW Hydrae (TW Hya). This high density in TW Hya was 
also observed using the XMM-Newton Reflection Grating Spectrometer (Stelzer & Schmitt 
2004). High n e at relatively low T e has now been measured in a number of CTTS systems 
(Gudel & Naze 2010, and references therein). 

Current address: Instituto de Astronomia y Fisica del Espacio, (IAFE), Buenos Aires, 
Argentina 
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TW Hya is nearby (57 pc) and, with an inclination angle of 7° (Qi et al. 2004), 
provides a pole-on view of the star where accretion occurs. Chandra observed TW Hya 
with the HETG for 489 ks in 2007 to obtain a spectrum with high signal-to-noise ratio. 
The spectrum provides a number of diagnostic line ratios not yet obtained for other CTTS 
systems (Brickhouse et al. 2010, hereafter Paper I). In particular, diagnostics indicate 
three emission components: a hot corona (~10 MK), an accretion shock (2.5 MK) in good 
agreement with model predictions, and a large post-shock region (<2 MK) with 30 times 
the mass of the shock front. O VII line ratio diagnostics of this third component rule out 
its origin in the radiative cooling zone of the shock and suggest instead its production by 
accretion-driven heating of surrounding stellar atmosphere. 

In this Letter we present new results from the 2007 Chandra spectra of TW Hya, 
focusing on the accretion shock component. Line ratio diagnostics measured from three 
different pointings show significant variation. Ne IX line ratio diagnostics constrain the 
fundamental parameters of a basic accretion model. In particular, the absorbing column 
density Njj derived from Ne IX line ratios provides a new diagnostic of the accretion shock. 



2. Spectral Analysis 

Chandra HETG spectra of TW Hya are analyzed from three different pointings 
spanning 2007 February to 2007 March — observation identifications (obsids) 7435 
(153.3 ks), 7437 (157.0 ks), and 7436 (158.4 ks). A fourth observation of TW Hya (obsid 
7438) is too short to provide sufficient signal-to-noise ratio for our purposes here. 

Table 1 presents fluxes from the three pointings for the strongest observed lines of 
Mg XI, Ne X, Ne IX, O VIII, and Fe XVII. In Paper I we argued that these ions are formed 
in the cooling column of the shock, whereas lines from Mg XII, Si XIII, and Si XIV are 
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formed in a stellar corona and lines from O VII are produced in a large post-shock region 
associated with the soft excess found in other CTTS (Gudel & Telleschi 2007; Robrade & 
Schmitt 2007). The Ne X lines are likely to be formed in both the accretion shock and the 
hot corona. We include them here in order to assess this assumption. 

The He-like Ne IX system provides line ratio diagnostics for T e , n e , and N H (Paper I). 
Fortunately, the atomic data for this ion are now accurate for diagnostic studies (Chen et al. 
2006). Figure 1 shows the theoretical T e -dependent curve for Ne IX (Foster et al. 2012) with 
the measured G-ratios [(forbidden plus intercombination) /resonance] overplotted. During 
the third observation T e is significantly higher than during the first two observations. This 
change in T e provides one of the primary motivations for this study. 

The observed Ne IX series ratios rule out resonance scattering, consistent with the 
measured line broadening (Paper I). The line optical depth is reduced by two orders of 
magnitude due to the modeled filling factor (Table 2) and turbulent broadening. Figure 2 
shows Nh inferred from the measured Ne IX resonance Hea (A13.45) /He/3 (All. 54) 
ratios, using a neutral to near-neutral photoelectric absorber. N H is smaller in the third 
observation compared with the first two, providing additional motivation for this study. To 
derive Nh from the Hea/He/3 line ratio, the Independence is first determined from the 
observed G-ratio. Theoretical curves of Hea/He/3 are then plotted as functions of Nh, for 
different values of T e ranging from 1.6 to 4.0 MK. Table 2 gives the derived values for T e , 
N H , and n e . We note that n e derived from the R-ratio (forbidden/intercombination; not 
shown) is only marginally significantly different in the third pointing. Since the values of 
n e from all three pointings are more than an order of magnitude higher than observed from 
purely coronal sources at Ne IX temperatures (Sanz-Forcada et al. 2003), the assumption 
that Ne IX forms primarily in the accretion shock remains justified. 

Fe XVII also originates primarily in the shock (Paper I). We use the observed A15.01 
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("3C") to A15.26 ("3D") line ratios to obtain T e « 3.0, 2.1, and 4.6 MK for pointings 1, 2, 
and 3, respectively. The temperature sensitivity of 3C/3D arises from an inner shell Fe XVI 
line blend with 3D (Brown et al. 2001). While these T e values are not entirely consistent 
with those from Ne IX, they support the trend for higher T e during the third pointing. 
Notably, the ratio for the third observation is consistent with the absence of Fe XVI, as 
expected at the higher temperature. 

We can test the overall consistency of the derived parameters by comparing the 
observed line fluxes with predictions. For a single component model, the observed 
line intensity I obs = W\ A e(T e ) EM/(AttR 2 ), where W\ is the wavelength-dependent 
transmission through the absorber (Morrison & McCammon 1983; Wilms et al. 2000), A 
is the relative elemental abundance, e is the emissivity from AtomDB (Smith et al. 2001; 
Foster et al. 2012), EM is the emission measure {EM = n 2 e V ', where V is the volume) and 
R is the distance to the source. The absorption model assumes solar abundances, which 
could introduce an error as large as a factor of two in Nh- 

We use the summed observation to determine A and then predict fluxes for the three 
pointings. Abundance ratios for oxygen, magnesium, and iron compared to neon are in 
the range of relative abundances found in Paper I. The observed fluxes of Ne IX A13.45 
along with the values of n e , T e , and W\ determined from the Ne IX diagnostics, give V 
for each pointing. The predicted fluxes of the resonance lines Mg XI A9.17, O VIII A18.97 
and Fe XVII A15.01 are reasonably consistent with observations; however, the observed 
fluxes show smaller differences than predicted. For all three ions, a smaller difference in the 
transmission W\ gives better agreement. For example, decreasing N H from ~ 3.0 x 10 21 
cm" 2 to ~ 2.2 x 10 21 cm" 2 brings the O VIII line for the first two pointings into good 
agreement, while increasing Nh from 0.9 xlO 21 cm" 2 to a value between 1.0 and 1.3 xlO 21 
cm -2 brings the line fluxes for all three ions into excellent agreement for the third pointing. 
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The consistency for Mg XI and Fe XVII is somewhat worse for the first two pointings, 
probably because their emissivity functions are steeply dropping at such low T e . Given 
order of magnitude differences in emissivities and transmission over the range in T e and 
Nh, the good agreement among the different shock diagnostics is remarkable. 

We predict Ne X Lya in a similar manner. Our predictions are low by a factor of ten 
for the first two pointings, but by only a factor of two for the third. Ne X arises primarily 
from the corona, but has a larger accretion contribution during the third pointing compared 
with the other two. 



3. Determining the Parameters of the Accretion Shock 

The observed variation of Nh provides new support for the idea that the accreting 
pre-shock stream is the absorber of the shocked, X-ray emission (e.g. Lamzin 1999; 
Gregory et al. 2007; Paper I), since other potential sources of absorption (e.g. interstellar, 
circumstellar) are not expected to vary significantly. Moreover, the nearly face-on view of 
the TW Hya system makes the disk itself an implausible absorber. To explore what the 
pre-shock gas as absorber implies about the geometry of the system, we use the strong 
shock condition to obtain the electron density of the pre-shock gas (n e /4) from the observed 

irough the same 



density n e of the shock. We assume the observed column density Nh is t 
pre-shock gas, and estimate the path length < I >= AN H /n e . From Ne IxQ N H = 1.8 x 10 21 
cm -2 and n e = 3.2 x 10 12 cm -3 , and thus < I >= 2.2 x 10 4 km. This path length is much 
smaller than the stellar radius, and thus geometrical effects due to the curvature of the 
accretion column are insignificant. This estimate of < I > is also an order of magnitude 



2 These values are from the summed spectrum, with a small correction for the In- 
dependence of the Ne IX R-ratio, following Smith et al. (2009). 
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larger than the length of the post-shock cooling zone (derived in Paper I). Absorption by 
the cooling column appears to be ruled out. Absorption by intervening stellar atmosphere 
cannot be ruled out, although recent hydrodynamic models produce an "observable" shock 
in the upper chromosphere for conditions similar to those of TW Hya (Sacco et al. 2008; 
2010). Thus we proceed with the assumption that the pre-shock stream is the only absorber. 

In the standard model of magnetospheric accretion (Calvet & Gullbring 1998; Gunther 
et al. 2007), the shock temperature can only vary if the accreting material originates at 
different distances from the star. Konigl (1991) gives an expression for the inner truncation 
radius of the disk that depends only on the mass accretion rate M acc and the surface 
poloidal magnetic field strength B*, in addition to the stellar mass M* and radius -R*. Thus 
the observed changes in T e indicate variations in M acc , B*, or both. 

Together T e , n e , and N H constrain the fundamental parameters of the standard 
accretion model, namely M acc , B*, and filling factor /. We assume a solar- type magnetic 
dipole, aligned with the rotation axis of the star (Cranmer 2008, 2009). Material originates 
in the gaseous accretion disk between the inner truncation radius r in and an outer radius 
r out , both of which are dependent variables in the model. The accretion stream shocks near 
the surface of the star, forming a circumstellar ring at a latitude determined by the dipole 
field geometry and r in . We have run a cube of 50 x 50 x 50 models, independently varying 
-B*, and / for comparison with observations. The value of M acc is calculated for the 
whole star, i. e. in the pole-on case the model accretion rate is twice what we observe. The 
model parameters are given in Table 2 and illustrated in Figure 3 for the case where the 
free fall velocity is calculated from r in . Similar results are found if we compute the free fall 
velocity as the square root of the product of the velocities derived from material released at 
the inner and outer edges of the disk, respectively. 

We illustrate the model constraints using the diagnostic lines of Ne IX, since Ne IX is 
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the only He-like ion that produces lines strong enough to determine T e , n e , and Nh from 
the accretion shock in all three pointings. For the summed data and for each pointing, 
approximately 100 models, or 0.1% of the 125,000 models computed, are found to agree to 
within la of the derived T e , n e , and Nh- The summed spectrum, gives the average values 
M acc = 1.5 x 10~ 9 M Q yr _1 , = 600 G, and / = 3.5%. The accreting material originates 
at < 31?*. 

Here we have assumed a single T e and n e model for the shock. In Paper I we found / 
= 1.1% and M acc = 4 x 10~ 10 M yr _1 by computing T e , n e and the resulting line emission 
as functions of depth in the post-shock cooling zone, and assuming a truncation radius of 
4.5 -R*. We then integrated the emission over the cooling zone and compared to diagnostic 
line fluxes from the summed spectrum. The shock front temperature was 3.4 MK, compared 
with 2.5 MK here, while the shock-front electron density was slightly smaller at 2.3 x 10 12 
cm" 3 . The higher accretion rate derived here arises primarily from the three times larger /, 
with somewhat larger n e offset by somewhat smaller T e . 

We compare the derived M acc = 1.5 x 10~ 9 M yr" 1 to 0.4, 1.3, and 5 x 10~ 9 M yr" 1 , 
from Muzerolle et al. (2000), Donati et al. (2011), and Batalha et al. (2002), respectively. 
Curran et al. (2011) report consistent rates from optical diagnostics. Our value agrees 
well with these estimates. Donati et al. (2011) demonstrate using spectropolarimetric 
signatures that a strong quadrupole magnetic field (2.5 to 2.8 kG) dominates the magnetic 
structure in the photosphere of TW Hya; they also report a weaker dipole component (400 
to 700 kG) which is capable of disrupting the accretion disk at 3 to 4 R*, also consistent 
with our results. Several authors report significantly lower mass accretion rates from X-ray 
observations: 1.8 x 10~ 10 (Curran et al. 2011); 2 x HT 11 (Stelzer & Schmitt 2004); 2 x 10~ 10 
(Giinther et al. 2007); and, 6 to 7 x 10~ n (Giinther 2011) M yr" 1 . All of these studies 
use a lower N H value determined from global fits to the spectra, accounting for a factor of 
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three or larger difference between their accretion rates and ours. The remaining differences 
may be due to differences in absolute abundances, filling factors, or relative accretion and 
coronal contributions to the spectrum, and to different spectra and methods. 

The accretion parameters derived from the third pointing are strikingly different from 
those derived from the first two pointings. M acc drops by a factor of 5 while / is smaller by 
a factor of 7. The truncation radius moves farther out by about 70%. While the derived 
r in values may be underestimated by the single component shock, r in cannot be too large 
or the shock temperature would not show any measurable variation. Accounting for the 
systematic bias introduced by the single component simplification suggests an average r in 
value of 4 i?* from the summed spectrum, with a range from 3.4 to 6.1 i?* for the three 
pointings. 

The large variations in M acc and / are unlikely to be the result of rotation or viewing 
angle, and thus it appears that the accretion rate changes substantially while the stellar 
magnetic field changes only slightly, if at all. The single spectrum obtained with the MIKE 
double-echelle spectrograph on Magellan during the third Chandra pointing (Dupree et al. 
2012) shows a decrease in the average blue veiling to r = 0.8 from earlier values of r > 1.2 
during the second Chandra pointing, consistent with a decrease in the accretion rate. The 
values of M acc from the three pointings, 2.8, 2.2, and 0.53 x 10~ 9 M Q yr" 1 , are all well 
within the literature values, suggesting that at least some of the differences in the literature 
are due to intrinsic variability rather than to differences in measurement method. The 
observations indicate that at lower accretion rates, the stellar magnetic field exerts stronger 
control of the gas, producing a smaller hot spot at higher latitude, as expected. 

Episodic events at lower latitudes occur in the data of Donati et al. (2011), but on 
time-scales much shorter than the time-resolution of the X-ray spectroscopy. Veiling and 
spectral line diagnostics for both accretion and wind showed variability during ground-based 
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observations contemporaneous with the Chandra observations (Dupree et al. 2012). The 
time resolution of these variation signatures is of the order of several minutes, again much 
shorter than the effective time resolution of the X-ray spectroscopy. X-ray light curves 
derived from accretion lines show variability on timescales of a few ks (see also Paper I), 
although these variations are not correlated with the optical photometry, and an enhanced 
accretion event (during pointing 2) is followed by systematic changes in the optical emission 
lines (Dupree et al. 2012). A decrease in the accretion rate leads to a complex response by 
the X-ray light curve: a decrease in intrinsic emission measure is countered by an increase 
in observed flux due to a decrease in the soft X-ray absorption by the accreting column 
itself. 

4. Conclusions 

Significant variability in X-ray accretion diagnostics occurs during the Chandra HETG 
observations of TW Hya. In particular, variability of the column density N H implicates 
the pre-shock accreting gas as the X-ray absorber. Nh and n e measurements from He-like 
Ne IX lines support an absorbing path length short compared to the stellar radius but 
long compared to the shock cooling depth. With T e , n e , and N H together the X-ray data 
constrain all the fundamental parameters of the standard accretion model, namely M acc , 

and /. 

The resulting average mass accretion rate during the Chandra observation, 
M acc = 1.5 x 10~ 9 M yr _1 , based on Ne IX diagnostics alone, agrees well with optical and 
ultraviolet techniques. The observations indicate that the mass accretion rate decreases by 
a factor of 5 during one pointing of ~150 ks, compared with two others. This accretion rate 
variation appears to arise from changes in the inner radius of the accretion disk. 
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Table 1. Emission Line Fluxes 



Ion A re / a Flux (uncorrected for absorption) 
(A) (l(r 6 ph cm- 2 s- 1 ) 



Pointing l b Pointing 2 b Pointing 3 b 



Mg XI 


9.17 


2.5 


± 


0.4 


1.7 


± 


0.5 


2.4 


± 


0.4 


NeX 


9.71 


1.8 


± 


0.8 


3.1 


± 


0.5 


2.8 


± 


0.6 


NeX 


10.24 


7.0 


± 


0.7 


9.0 


± 


2.0 


8.7 


± 


0.8 


i\Tq ty 
iNe ia. 


1 1 nn 
11. UU 


( .6 


1 

± 


U. ( 


0.0 


1 

± 


l.U 




1 

± 


l.U 


Ne IX 


11.54 


19.9 


± 


1.5 


23.7 


± 


1.8 


24.2 


± 


1.6 


NeX 


12.13 


63.8 


± 


3.1 


74.1 


± 


3.3 


77.2 


± 


3.4 


Fe XVII 


12.27 


3.2 


± 


0.8 


4.7 


± 


1.0 


2.1 


± 


0.6 


Nc IX 


13.45 


128.3 


± 


6.9 


144.7 


± 


7.0 


187.8 


± 


8.0 


Ne IX 


13.55 


93.5 


± 


5.5 


97.8 


± 


5.9 


115.4 


± 


6.2 


Ne IX 


13.70 


47.4 


± 


4.0 


50.4 


± 


3.4 


53.9 


± 


4.3 


Fc XVII 


15.01 


28.6 


± 


3.5 


29.2 


± 


3.5 


43.4 


± 


4.1 


Fc XVII 


15.26 


11.7 


± 


2.5 


15.0 


± 


3.2 


15.8 


± 


2.4 


VIII 


16.01 


27.1 


± 


3.2 


24.3 


± 


3.8 


32.7 


± 


3.9 


Fe XVII 


16.78 


25.6 


± 


4.7 


22.4 


± 


4.0 


25.3 


± 


4.3 


Fe XVII 


17.05 


25.5 


± 


7.6 


28.4 


± 


4.5 


31.8 


± 


4.3 


Fe XVII 


17.10 


23.0 


± 


6.4 


28.7 


± 


4.6 


29.4 


± 


4.2 


VIII 


18.97 


201.3 


± 


15. 


210.1 


± 


15. 


201.0 


± 


16. 



"■Reference wavelengths from AtomDB v2.0 (Foster et al. 
2012). Multiplets are intensity- weighted averages. 

b Pointings 1, 2, and 3 correspond to mid times 48.0, 58.5, 
and 61.2 (JD-2454100), respectively. 
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Table 2. Physical Conditions and Accretion Model Parameters 



Parameter Pointing l a Pointing 2 a Pointing 3 a Total b 

T e (MK) C 1.9±g;| 2.3^1 3.lt°;J 2.5±g| 

n e (10 12 cm- 3 ) c 3.0tg;£ 3.1±8:| 3.9±g"i 3.2t°; 2 

7V ff (10 21 cm- 2 ) d 3.2tJ-i 2.8±g:l 0.9t°j 1.8±g;« 

M acc (10~ 9 M yr- 1 ) 2.78 (2.02-3.56) 2.19 (1.68-2.95) 0.53 (0.19-1.15) 1.45 (1.05-1.84) 

Poloidal B* (G)° 523. (389.-679.) 628. (456.-797.) 770. (359.-1510.) 614. (494.-797.) 

Filling Factor (%) c 8.19 (5.08-12.3) 5.69 (3.91-8.63) 0.95 (0.29-2.46) 3.46 (2.52-4.59) 

r m /R, 1.77 (1.58-1.99) 2.11 (1.84-2.42) 3.60 (2.58-4.89) 2.34 (2.07 2.60) 

w/R, c 2.23 (1.86-2.67) 2.57 (2.10-3.04) 3.82 (2.67-5.33) 2.68 (2.33-3.14) 

Latitude f 44.6° 49.0° 58.7° 50.8° 

a Pointings 1, 2, and 3 correspond to mid times 48.0, 58.5, and 61.2 (JD-2454100), respectively. 

b Derived properties T e , n e and Nh are from Brickhouse et al. (2010), except for a minor revision 
to n e , for which we have taken the Independence of the n e diagnostic into account here. The error 
on Nh is estimated in this work. 

c Errors on T e and n e use la errors on the observed line ratios. 

d Errors on Nh include la errors on T e (from the G-ratio) as well as la errors on the Hea/He/3 ratio. 

Values given are the average values from the models that meet the criteria for each time segment, 
described in the text. The range in parentheses gives the minimum and maximum values for the same 
set of models. 



f Average latitude of the field line which the accretion stream follows. 
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1.5xl0 6 2.0xl0 6 2.5xl0 6 3.0xl0 6 3.5xl0 6 4.0xl0 6 

T e (K) 

Fig. 1. — Theoretical G-ratio (i.e., ratio of forbidden plus intercombination to resonance line 
fluxes) for Ne IX as a function of T e (solid curve) from AtomDB (Foster et al. 2012; Chen 
et al. 2006). Overplotted are the observed ratios with lcr errors (vertical error bars) for the 
individual observations (open circles, numbered in time order) and for the summed spectrum 
(solid circle). Horizontal error bars (dashed lines) show how the errors on the observed ratios 
translate to the errors on T e . 
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Fig. 2. — Ratio of the flux of the resonance line (Hea) to that of He/3 as a function of 
hydrogen column density (dashed curves), for given log [T e ] (K) from 6.2 to 6.6 as labeled. 
Absorption models are from Morrison & McCammon (1983) and line emissivities are from 
Foster et al. (2012). Overplotted are the observered ratios with la errors (vertical error bars) 
for the individual observations (open circles, numbered in time order) and for the summed 
spectrum (solid circle). Observed points are placed on the curve appropriate for their T e 
as determined from their G-ratios. Horizontal error bars (dotted lines) show how the errors 
both on the observed Hea/He/? ratios and on T e (determined from the G-ratios) translate 
to the errors on N H . 
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Fig. 3. — Parameter ranges for the dipole accretion model, constrained by Ne IX diagnostics 
from time segments 1 (JD 48.0; red crosses), 2 (JD 58.5; green triangles), 3 (JD 61.2; blue 
diamonds), and the summed data (black X's). The model parameters plotted vs the inner 
disk radius are outer disk radius (a), mass accretion rate (b), filling factor at the stellar 
surface (c), and polar magnetic field at the star (d). 



